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ABSTRACT 

We develop a hybrid galaxy formation model which uses outputs from an N-body 
simulation to follow the merger histories (or "merger trees" ) of dark matter halos and 
treats baryonic processes, such as the cooling of gas within halos and subsequent star 
formation, using the semi-analytic model of Cole et al. We compare this hybrid model 
to an otherwise identical model which utilises merger tree realisations generated by a 
Monte-Carlo algorithm and find that, apart from the limited mass resolution imposed 
by the N-body particle mass, the only significant differences between the models are 
due to the known discrepancy between the distribution of halo progenitor masses 
predicted by extended Press Schechter theory and that found in N-body simulations. 
We investigate the effect of limited mass resolution on the hybrid model by comparing 
to a purely semi-analytic model with greatly improved mass resolution. We find that 
the mass resolution of the simulation we use, which has a particle mass of 1.4 x 
10^*^/1"^ Mq, is insufficient to produce a reasonable luminosity function for galaxies 
with magnitudes in the bj band fainter than -17. 
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1 INTRODUCTION 

Hierarchical models of galaxy formation must describe both 
the growth and collapse of density perturbations to form 
dark matter halos and the baryonic processes which lead to 
the formation of stars. Despite uncertainty as to the exact 
nature of the dark matter itself, the formation and evolution 
of dark matter halos appears to be reasonably well under- 
stood. The two main approaches to this problem are direct 
numerical simulations and analytic techniques such as the 
Press-Schechter theory (Press & Schechter 1974). Encour- 
agingly, the mass functions of dark matter halos predicted 
using these very different approaches are found to agree to 
within 50% (Gross et al. 1998; Governato et al. 1999; Jenk- 
ins et al. 2001). The analytic model described by Sheth & 
Tormen (2001) based on the assumption that objects col- 
lapse ellipsoidally rather than spherically achieves even bet- 
ter agreement with N-body simulations. Mo & White (2002) 
present halo abundances from this and several other models. 

This understanding of the hierarchical build up of struc- 
ture provides the starting point for semi-analytic models of 
galaxy formation, which attempt to follow the development 
of galaxies from primordial density fluctuations. In semi- 
analytic models, merger histories for dark matter halos may 



be taken directly from dark matter simulations (e.g. Kauff- 
mann et al. 1999; van Kampen et al. 1999). Alternatively, 
extensions to the Press-Schechter theory which predict the 
conditional halo mass function (Bond et al. 1991; Bower 
1991) and halo survival times, formation times and merger 
rates (Lacey & Cole 1993) may be used to construct reali- 
sations of merger histories for individual halos. Simple an- 
alytic modelling is then used to follow the evolution of the 
baryonic component, including prescriptions for processes 
such as star formation and its possible effects on the re- 
maining gas. Semi-analytic models (e.g. Cole 1991; Lacey & 
Silk 1991; White & Frenk 1991; Cole et al. 1994; Somerville 
& Primack 1999; Cole et al. 2000) have successfully repro- 
duced many observable properties of galaxies, such as the 
local field galaxy luminosity function and distributions of 
colour and morphology. When combined with N-body simu- 
lations, semi-analytic models have also successfully repro- 
duced galaxy clustering properties (e.g. Governato et al. 
1998; Kauffmann et al. 1999; Benson et al. 2000; Wechsler 
et al. 2001). 

Semi-analytic models utilising merger trees generated 
using algorithms based on the extended Press-Schechter 
(EPS) formalism have two closely related advantages over 
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models which take merger histories from N-body simula- 
tions. Creating Monte-Carlo realisations of merger trees for a 
set of halos generally requires fewer computing resorces than 
carrying out an N-body simulation of a similar number of 
halos. In both cases, improving the mass resolution increases 
the computational load, but since the load is much less in the 
Monte-Carlo case, significantly better mass resolution may 
be achieved. Methods based on the Press-Schechter theory, 
however are only applicable to initially Gaussian fluctation 
fields. N-body simulations, on the other hand, have the ad- 
vantage that the non-linear evolution of density fiuctuations 
is followed in complete generality, without the need for any 
of the assumptions involved in creating EPS merger trees. 

There are advantages to both of these methods, and 
which is more appropriate depends on the problem being 
addressed. In this paper we investigate the effects of the 
choice of merger trees on the predictions of one particular 
semi-analytic model. We describe a new method of extract- 
ing merger trees from an N-body simulation and incorpo- 
rate these merger trees into a semi-analytic galaxy formar 
tion model based on that of Cole et al. (2000) . We compare 
the predictions of this model to those of a similar model 
utilising Monte-Carlo realisations of halo merging histories. 
In order to identify the reasons for the discrepancies that we 
find, we determine the changes that must be made to the 
Monte-Carlo model to reproduce the N-body results. 

The use of N-body merger trees in semi-analytic models 
allows a halo-by-halo comparison between the semi-analytic 
treatment of baryonic processes, such as gas cooling, and 
direct numerical simulations of galaxy formation. In a com- 
panion paper (Helly et al. 2002) we carry out such a compar- 
ison between a "stripped down" version of the semi-analytic 
model described in this paper and a smoothed particle hy- 
drodynamics simulation of a cosmological volume. 

This paper is laid out as follows. In Section 2 we explain 
how we obtain merger trees from an N-body simulation. 
In Section 3 we investigate the effect on our semi-analytic 
model of utilising merger trees derived from N-body simula- 
tions rather than Monte-Carlo realisations. In Section 4 we 
present our conclusions. 



2 EXTRACTING MERGER TREES 

We now present the method we used to calculate the merger 
histories of dark matter halos identified in an N-body simu- 
lation. The simulation, which will be referred to as the GIF 
simulation, was carried out by the Virgo Consortium using 
a parallel adaptive particle-particle/particle-mesh (AP'^M) 
code known as Hydra (Couchman, Thomas & Pearce 1995; 
Pearce & Couchman 1997) as part of the GIF project. 
The simulation assumes the ACDM cosmology with mean 
mass density parameter f2o = 0.3, cosmological constant 
Ao = 0.7 in units of power spectrum shape param- 

eter r = 0.21, present day rms linear fluctuation amplitude 
in 8/!.~^Mpc spheres erg — 0.90, and Hubble constant h = 0.7 
in units of 100kms~^Mpc~^. It contains 256^ dark matter 
particles each of mass 1.4 x 10^°/i"^Mo in a box of side 
141.3/i^^Mpc. The gravitational softening length in the sim- 
ulation is 30/i~^kpc at « = 0. This simulation is described 
in more detail by Jenkins et al. (1998), where it is referred 
to as ACDM2, and by Kauffmann et al. (1999). While halo 



catalogues and merger trees based on this simulation are 
publically available, here we make use of only the simula- 
tion outputs themselves and construct merger trees using a 
somewhat different algorithm to that of Kauffmann et al. We 

use 44 output times from the simulation which are spaced 
equally in log^Q(l -|- z) between z = Q and z ^ 20. 

2.1 Identifying Halos 

In order to construct merger histories for dark matter halos 
in an N-body simulation, a catalogue of halos must be pro- 
duced for each simulation output using a group finding al- 
gorithm. The algorithm used here is the "friends of friends" 
(FOF) method of Davis et al. (1985), which simply links to- 
gether any particles with separations less than the linking 
length b, usually expressed in terms of the mean interparti- 
cle separation. Given sufficiently large numbers of particles 
in each object, the FOF algorithm finds regions bounded 
by a surface of constant density. The density threshold is 
proportional to 1/6^. 

The FOF approach has the advantage that it imposes 
no constraints on the geometry of the halos identified, but it 
may occasionally artificially join two nearby halos if a tran- 
sient "bridge" of a few particles forms between them. It will 
be seen in Section 2.2 that this can cause problems when 
attempting to generate merger trees using FOF group cata- 
logues, and a method of identifying and splitting artificially 
joined halos is described in Section 2.2. 

The usual choice for the linking length in cosmologies 
with n = 1 is 6 = 0.2 (e.g. Lacey & Cole 1994), which 
identifies halos with a mean density similar to that predicted 
by the top hat spherical collapse model (Cole & Lacey 1996). 
However, in cosmologies with n < 1 there is no rigorous 
justification for any particular choice. Here, we choose to 
set h = 0.2 at all redshifts as in the = 1 case. See Eke, 
Cole & Frenk (1996) and Jenkins et al. (2001) for further 
discussion. 

The other parameter needed by the FOF algorithm is 
the minimum number of particles, A^min, required to consti- 
tute a group. It is important that Nmin be as small as pos- 
sible, since detailed merger trees can only be obtained for 
halos much larger than the smallest resolvable group. Kauff- 
mann et al. (1999) found that in their simulations groups as 
small as 10 particles are dynamically stable systems and that 
for 95% of these groups, 80% of the particles remain in the 
same group at subsequent times. 

We therefore identify halos using a linking length h = 
0.2 at all redshifts, with a minimum group size of ten parti- 
cles. The resulting catalogues may still contain some groups 
which consist of unbound particles which happen to be close 
together at this particular timestep. To remove these, we fol- 
low Benson et al. (2001b) and calculate the total energy of 
each group. Unbound groups are not immediately discarded, 
because they may only be unbound due to the presence of a 
small number of fast moving particles. The binding energy 
of each particle is calculated, and the least bound particle 
is removed from the group. This is repeated until the group 
becomes bound. If half of the particles are removed or the 
group is reduced to less than Amin particles we discard it. 
Up to 5% of all groups are discarded, with a similar num- 
ber of groups being reduced in mass by this procedure. The 
affected groups generally consist of around 10-20 particles. 
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We use the procedure described above to generate halo 
catalogues for 44 simulation outputs between redshifts z = 
20 and 2 = 0, spaced approximately evenly in logj()(l + z). 

2.2 Constructing N-body Merger Trees 

In an idealised picture of the process of hierarchical structure 
formation (e.g. Press-Schechter theory), dark matter halos 
may increase in mass by mergers, but cannot lose mass. 
Consequently, any halo identified in a simulation prior to 
the final output time should still exist at subsequent out- 
put times, although it may have become subsumed within 
a larger halo through a merger. In any case, the constituent 
particles of the original halo should still all be members of a 
single group. It should therefore be possible to identify each 
halo in the simulation as a progenitor of a single halo at the 
next output time. 

In practice there are several ways in which a halo can 
lose particles. Halos may be disrupted by tidal forces caused 
by other nearby halos. The masses of simulated halos can 
also fluctuate because the FOF algorithm imposes a some- 
what arbitrary boundary on the halo and outlying particles 
which are considered group members at one timestep may 
lie just beyond the boundary at the next timestep. 

The technique we use to determine merger histories is 
intended to take into account this uncertainty in the def- 
inition of a halo and a possible loss of particles. First, we 
consider two adjacent output times from the simulation, ti 
and t2, where t\ < t2- Each halo at time ti is labelled as a 
progenitor of whichever halo at time t2 contains the largest 
fraction of its particles. This process is repeated for all pairs 
of adjacent output times. It is then straightforward to trace 
the merger history of each halo which exists at the final out- 
put time. Fig. 1 shows an example of a merger tree created in 
this way for a halo with a final mass of about 9 x 10^^ M© , 
or around 700 particles. 

In the semi-analytic model used here, galaxies are as- 
sumed to form at the centres of dark matter halos, so the 
centre of each halo in the merger tree must be defined. We 
choose to follow Kauffmann et al. (1999), who identified 
the most bound dark matter particle as the position of any 
galaxy which forms in the halo. We define the binding en- 
ergy of a particle as the sum of its kinetic energy and the 
gravitational potential energy due to the other particles in 
the halo. This approach differs from that of Benson et al. 
(2001a), who associated the central galaxy in a halo with the 
centre of mass. Once a galaxy forms it is assumed to follow 
this particle until the parent halo merges with another halo 
and dynamical friction, calculated as described in Cole et al. 
(2000), causes the galaxy to merge with the central galaxy 
of the new halo. We therefore check that the most bound 
particle of a halo remains a member of the same halo as the 
majority of the halo's constituent particles at the next out- 
put time. If this is not so, wo choose the most bound particle 
from those which are in the correct halo at the later output 
time. This problem generally only occurs in smaller halos 
which may be easily disrupted. 

During the construction of the merger trees, we also at- 
tempt to deal with the problem mentioned in Section 2.1 — 
the possibility that nearby halos may be artificially linked 
by the FOF algorithm. The problem occurs if two halos be- 
come temporarily linked by a transient "bridge" of parti- 
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Figure 1. An example of a merger tree obtained from the GIF 
simulation for a halo of mass 9 X 1012/i-1Mq at rcdshift z = 0. 
Each circle represents a dark matter halo identified in the simu- 
lation, the area of the circle being proportional to the halo mass. 
The vertical position of each halo on the plot is determined by 
l°Sio(l + the redshift at which it exists, the horizontal po- 
sitioning is arbitrary. The solid lines connect halos to their pro- 
genitors. The solid line in the panel on the left-hand side shows 
the fraction of the final mass contained in resolved progenitors 
as a function of redshift. The dotted line shows the fra<;tion of 
the final mass contained in the largest progenitor as a function of 
redshift. 



cles which causes the FOF group finder to consider them 
as a single, large group. When the bridge is later broken, 
the group splits, leaving the two original halos. Our tree 
building method would identify the large, joined group as a 
progenitor of the larger of the two final groups. 

These situations are identified by looking for groups at 
the earlier time ti whose particles are shared between two or 
more groups at the subsequent output time t2. This indicates 
that between times ti and t2 the group has split into smaller 
groups which we refer to here as "fragments" . 

We split such spuriously joined groups into one new 
group for each fragment which contains more than A^min of 
its constituent particles. Particles belonging to one of these 
fragments at time t2 are assigned to the corresponding new 
group at the earlier time ti. Particles belonging to no frag- 
ment, or to a fragment with fewer than A'min particles from 
the joined group, are assigned to the new group correspond- 
ing to the fragment "closest" to their position at time ti. 
The separations used are weighted by a factor M~^^^ to ac- 
count for the spatial extent of the groups, where M is the 
mass of the fragment. 

The splitting procedure is first carried out for halos at 
the penultimate timestep and then repeated for each earlier 
output time in order of increasing redshift. For each timestep 
a modified group catalogue is produced, which is then used 
to determine whether any halos at the previous timestep 
need to be split. This ensures that if any bridge between a 
pair of halos persists for more than one timestep the halos 
are split at each timestep where the bridge exists. 
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2.3 Mass Conservation 

In the GALFORM somi-analytic model of Cole et al. (2000), 
halos may gain mass through mergers with other halos. The 
mass of a halo always increases with time, and the difference 
between the mass of a halo and the sum of the masses of its 
progenitors is due to the accretion of small, unresolved dark 
matter halos. 

The N-body merger trees may contain halos which de- 
crease in mass from one timestep to the next for the reasons 
described in Section 2.2 — the nature of the definition of 
a halo imposed by the FOF group finder and the possibil- 
ity of disruption by tidal forces. Consequently, a halo in a 
N-body merger tree may be somewhat less massive than 
its progenitors. In the GALFORM model this corresponds to 
the unphysical situation where a negative amount of mass 
is accreted in the form of sub-resolution halos. 

The solid lines in Fig. 2 show the distribution of the 
ratio EMprog/Mhaio, where Mhaio is the mass of a halo and 
EMprog is the total mass of the immediate progenitors of 
the halo, which exist at the previous timestep. Halos at all 
timesteps (other than the first) are included. If these merger 
trees had been created using the technique of Cole et al. 
(2000), then this ratio would always be less than one. It 
can be seen from Fig. 2 that for halos less massive than 
about 10^^/i~^Mo the total mass in the progenitors can oc- 
casionally exceed the mass of the halo they form at the next 
timestep by up to 50%. More massive halos are less affected, 
but there are still rare instances where the largest halos have 
progenitors with masses 5-10% greater than the mass of the 
halo. 

Mass conservation can be forced on the N-body merger 
trees by simply adjusting the masses of some of the halos. 
Two opposite approaches to the problem are possible. Mass 
can be added to those halos which axe less massive than their 
progenitors, or mass can be removed from the progenitors 
themselves. In order to show that the changes made to the 
halo masses have little effect on the semi-analytic model, we 
create merger trees using both methods. 

Enforcing the conservation of mass in merger trees by 
adding mass is relatively straightforward. If a halo is less 
massive than its progenitors, its mass is increased to match 
that of the progenitors. The halo may, in turn, be a progeni- 
tor of a later halo which may now become less massive than 
its own progenitors. This later halo's mass will then also be 
increased. Changes made to halo masses at early times may 
therefore propagate to later times. 

Similarly, if mass is removed from a halo to force con- 
servation of mass, it may become less massive than its pro- 
genitors and reductions in mass could then propagate to 
earlier times. We attempt to remove mass in such a way 
as to minimize the effects on earlier halos. Each halo has 
a certain amount of "excess" mass beyond that of its pro- 
genitors, which was accreted over the last timestep in the 
form of sub-resolution objects. This mass, if it exists, may 
be removed without the change propagating to earlier ha- 
los. When a halo which is less massive than its progenitors 
is found, mass is first removed from the excess mass of the 
largest progenitor. If still more mass must be removed, it 
is taken from the excess mass of the other progenitors in 
decreasing order of mass. If all of the excess mass of the 
progenitors is removed and yet more mass needs to be taken 



away, the masses of all of the progenitor halos are simply 
scaled down by a constant factor. 

The dotted lines in Fig. 2 show the sizes of the changes 
we are forced to make when we enforce mass conservation 
by adding mass to halos. These lines show the distribution 
of the ratio EMprog/Mhaio if EMprog is evaluated after the 
progenitors of the halo at all previous timesteps have been 
made at least as massive as their own progenitors. Mhaio is 
still the original halo mass. Where this ratio exceeds 1, it is 
the factor by which Mhaio must be scaled to ensure that the 
halo is at least as massive as its progenitors. It can be seen 
that the required changes to individual halos are generally 
small, and adjustments are required much less frequently in 
well resolved halos. However, the masses of a minority of ha- 
los are affected quite significantly and it is necessary to show 
that these changes do not affect the galax;y population pre- 
dicted by the semi-analytic model. The algorithms described 
above are two opposite ways of dealing with the problem of 
mass conservation in the merger trees. While artificially al- 
tering the halo masses is clearly not ideal, if, as is the case, 
both methods produce very similar results when the merger 
trees are fed into the semi-analytic model we can then con- 
clude that the changes we have made are insignificant. This 
comparison is carried out in Section 3.3. 



3 COMPARISON BETWEEN GALFORM AND 
N-BODY GALFORM 

In this section we describe our semi-analytic model, indi- 
cating how it differs from the model of Cole et al. (2000) on 
which it is based. We also explain how merger trees obtained 
from a simulation may be incorporated into the model. 

3.1 The N-body GALFORM model 

We use the GALFORM semi-analytic model to treat the pro- 
cess of galaxy formation within the dark matter halos in the 
GIF simulation. The model is described in detail by Cole et 
al. (2000) so here we present only a brief description of fea- 
tures that arc important to this work. The original model of 
Cole et al. will be referred to as "standard galform" , and 
the version using merger trees taken from a simulation will 
be referred to as "N-body galform" . 

The starting point for the standard GALFORM model is a 
set of merger trees created using a Monte-Carlo technique. 
The history of each halo is divided into a number of dis- 
crete timesteps. Extended Press-Schechter theory is used to 
estimate the probability that a halo "fragments" into two 
progenitors when a step back in time of size St is taken. 
The masses of the fragments are chosen at random from a 
distribution consistent with extended Press-Schechter the- 
ory. Halos are repeatedly split in this way to create merger 
trees. A mass resolution limit is imposed on the merger trees, 
below which progenitors are considered to be material ac- 
quired through continuous accretion. The mass resolution is 
normally set sufficiently high that the results of interest are 
not sensitive to its value. In the N-body galform model, 
we replace these merger trees with those calculated directly 
from the GIF simulation as described in Section 2.2. The 
mass resolution limit is then determined by the mass of the 
smallest halo which can be resolved in the simulation. 
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Figure 2. The solid lines show the distribution of the ratio of the total mass of the immediate jHogcnitors of a halo, SMprog, to the 
mass of the halo at the next timestep, M^aio- Each panel shows the distribution of SMprog/Mimlo for halos in the mass range shown at 
the top of the panel. The dotted lines show the distribution of SMprog /Mhaio if SMprog is evaluated after the progenitors have been 
increased in mass to at least the total mass of their progenitors. Where this ratio is greater than 1, it is the factor by which Mhaio must 
be changed to ensure mass conservation if we choose to add mass. 



The dark matter halos in the merger tree are assumed 
to be spherically symmetric with the radial density profile 
of Navarro, Prenk & White (1996, 1997): 



r/rNFVf{r/rNFVf + 1)'' 

where tnfw is the scale radius of the halo and is related to 
the concentration parameter, c, defined by Navarro, Prenk 



& White (1997) through tnfw = Tviriai/c, where rviriai is 
the virial radius of the halo. The concentration parameter is 
set using the method described in the appendix of the same 
paper. We do not allow for any scatter in the concentration 
parameter as a function of halo mass. 

Our treatment of the cooling of gas within halos is iden- 
tical to that of Cole et al. (2000). Initially, the amount of 
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gas in each halo is taken to be equal to the mass of the halo 
times the universal baryon fraction. The gas is assumed to 
be shock-heated to the virial temperature of the halo when 
it forms. We assume that the radial density profile of the 
gas is given by 

Pgas(r) OC l/(r^ + reore), (2) 

where the core radius is given by rcoro/j'NFW ~ 1/3 in ac- 
cordance with the simulations of Navarro et al. (1995). This 
core radius is allowed to grow with time from an initial value, 
^TOrei as gas is removed by cooling in order to maintain the 
same gas density at the virial radius. This ensures that the 
pressure at the virial radius, which would be maintained by 
shocks from infalling material, remains unchanged. 

To determine the rate at which gas can cool and form a 
disk at the centre of the halo, the cooling time of the gas is 
calculated as a function of radius using the cooling function 
of Sutherland & Dopita (1993). Gas which has had time to 
cool and fall to the centre of the halo is added to the disk 
where it is available to form stars. 

When halos merge, the most massive galaxy becomes 
the central galaxy in the new halo. The resolution of the 
simulations used hero is insufficient to follow the evolution 
of substructure within the dark matter halos. Instead, the 
dynamical friction time scale, as defined by Lacey & Cole 
(1993), is used to determine when each satellite will merge 
on to the central galaxy. It should be noted at this point 
that the orbital parameters used to determine the dynam- 
ical friction time for each galaxy are assigned at random 
from a distribution consistent with the numerical results of 
Tormen (1997), even when using merger trees obtained from 
the simulation. 

3.2 Parameters in the N-body GALFORM model 

The GALFORM semi-analytic model requires a number of pa- 
rameters to be specified, which can be divided into three 
categories. There are numerical parameters, parameters de- 
scribing the background cosmology and parameters which 
describe the physical model of galaxy formation. 

The numerical parameters are the mass resolution, 
Mros, the number of timesteps in the merger tree and the 
starting redshift. In the N-body GALFORM model these are 
all constrained by the properties of the simulation used to 
obtain the merger trees. The mass resolution is the mass of 
the smallest halo which our group finding algorithm can re- 
solve, there is one timestep for each simulation output and 
the starting redshift is the redshift of the first output. The 
cosmological parameters fio, Ao, h, as, T and, in the case of 
a simulation with a baryonic component, Qb, axe also fixed 
by the simulation. 

The remaining parameters allow us to vary the treat- 
ment of the processes involved in galaxy formation. The pa- 
rameters we are interested in are: 

• '^'coro- the initial size of the core in the radial gas density 
profile, specified in terms of tnew (see eqn. 2). 

• The evolution of rcore with time. The radius rcore may 
be a fixed fraction of t-nfw or it may be allowed to increase 
with time as described in Section 3.1 

• /df : A factor by which the dynamical friction time scale 
for a satellite galaxy, which is used to determine when the 



galaxy merges with the central galaxy of the halo, may be 
scaled. Increasing /df reduces the rate at which galaxy merg- 
ers occur within halos. 

The other parameters in the model are the same as 
those in the reference model of Cole et al. (2000), with the 
following minor changes: tihot = 250kms~^ and /oiup = 0.5. 
The parameter Vhot determines the efficiency with which en- 
ergy injection from supernovae and young stars reheats and 
ejects cold gas from galactic disks. The parameter /eiup is 
used to decide the outcome of mergers between central and 
satellite galaxies. If the ratio of the mass of the satellite to 
the mass of the central galaxy is greater than /eiiip, any gas 
in the disks of the two galaxies is converted into stars and 
an elliptical galaxy is produced. If the ratio is smaller than 
/ciiip, any stars present in the satellite are added to the bulge 
of the central galaxy and any gas is added to the disk. These 
changes to the Cole et al. model are required to obtain a re- 
alistic luminosity function at a = with the higher baryon 
density, S7b = 0.038, which we use here. 

Our prescription for star formation differs slightly from 
that of Cole et al. In our model, the time scale for star 
formation is given by 

r. =ri'(I/disk/200kms-')"*, (3) 

where Vdisk is the circular velocity of the galaxy disk and the 
time scale, r°, is set to SGyr. We set a* = —2.5. The way 
Tt scales with redshift in this model results in reduced star 
formation and more gas rich mergers at high redshift and 
has been shown (Lacey 2002) to better reproduce the prop- 
erties of SCUBA and Lyman break galaxies. Kauffmann & 
Haehnelt (2000) also find that a star formation scheme with 
an increased star formation timescale at high redshift is re- 
quired to reproduce observations of damped Lya absorption 
systems and the increase in number density of bright quasars 
from X = to ^; = 2. It should also be noted that, for the 
purposes of this comparison, the details of our star forma- 
tion prescription are not critical, since the same scheme is 
used in both the standard and N-body galform models. 

3.3 Effects of mass conservation 

The upper panels of Fig. 3 show the galaxy luminosity func- 
tions in the bj and K bands predicted by the N-body gal- 
form model with the parameters of Section 3.2, using the 
two different methods described in Section 2.3 to enforce 
mass conservation in the merger trees. Over most of the lu- 
minosity range plotted, the two curves are essentially iden- 
tical but there appear to be more galaxies at very faint bj 
magnitudes when mass is removed from the merger trees. 
The majority of these galaxies formed in halos near the 
10 particle (~ 1.4 x lO^ft-^M©) mass resolution limit im- 
posed by the FOF group finder and their halos subsequently 
merged with other, larger dark matter halos. When mass 
conservation is enforced by removing mass from the merger 
trees (the dotted lines in Fig. 3) it is possible to end up with 
some halos with mass less than the resolution limit which 
can harbour galaxies with bj band magnitudes around -14 or 
fainter. If, instead, mass is added to halos less massive than 
their progenitors, then the merger trees contain no halos 
with masses below the FOF resolution threshold and hence 
fewer faint galaxies. 
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These sub-resolution halos often exist in the merger 
trees of larger halos and could affect the evolution of larger, 
brighter galaxies. However, the agreement of the luminosity 
functions suggests that any effect is insignificant. The global 
star formation history and TuUy-Fishcr relation shown in the 
lower panels of Fig. 3 are similarly unaffected. 

Overall, the choice of mass conservation method ap- 
pears to make very little difference to the quantities plotted 
in Fig. 3, which suggests that the small amounts of mass 
being added to or removed from the merger trees do not 
significantly affect the properties of the resulting galaxies. 
The only region of the luminosity function which is affected 
is largely populated by galaxies which formed in halos with 
little or no resolved merger history, where the model cannot 
be expected to give reliable results. For the remainder of 
this paper we choose to enforce mass conservation by adding 
mass to the merger trees since this does not introduce halos 
with masses below the resolution limit. 



3.4 Comparison with standard GALFORM 

The mass resolution of the merger trees taken from the GIF 
simulation is equal to 10 particle masses or 1.4x 10^^/i~^Mq, 
i.e. A/min = 10. This is much larger than the mass resolution, 
A/rcs = 5.0 X l(fh-^MQ, used by Cole et al. (2000). This 
will clearly affect the properties of the galaxies predicted by 
the N-body GALFORM model, since gas will be unable to cool 
and start forming stars until lower redshifts when halos with 
masses greater than Mres have formed. In order to investi- 
gate the effect of limited mass resolution on the N-body 
GALFORM model, we identify the properties of the merger 
trees which differ between standard and N-body GALFORM 
and use this knowledge to produce a modified version of the 
standard galform model which reproduces the behaviour 
of the N-body GALFORM model. We can then increase the 
mass resolution of the merger trees in the modified model 
and observe the effects on the predicted galaxy properties. 

There are four main reasons why the merger trees in 
the two models may differ. Firstly, there is the difference 
in mass resolution described above. Therefore, we initially 
degrade the mass resolution of the standard galform model 
to match that of the GIF simulation by setting the minimum 
halo mass, Mres , equal to the mass of (iVmin — 1) dark matter 
particles — any halo of this mass or less in the N-body 
simulation would not be identified by the EOF group finder 
and would not be included in the N-body merger trees. 

Secondly, Jenkins et al. (2001) have shown that the 
Press & Schechter (1974) halo mass function (used in the 
standard GALFORM model) differs somewhat from the mass 
function determined from N-body simulations. We replace 
the Press-Schechter mass function in the standard galform 
model with the mass function determined by Jenkins et al. 
This ensures that the distribution of halo masses at z = 
in the standard GALFORM model matches the distribution in 
the simulation. 

The number of timesteps also differs between the 
two models. In the standard galform model we use 150 
timesteps evenly spaced in logj^;, (1 + z), whereas in the N- 
body case we have only 44 simulation outputs. However, we 
find that if we degrade the time resolution of the standard 
GALFORM model to match that of the N-body model the 



properties of the galaxy populations predicted change very 
little. 

Finally, the distribution of progenitor masses for halos 
of a given mass predicted by the standard GALFORM model 
does not reproduce the distribution found in N-body simu- 
lations with complete accuracy. Benson et al. (2001a) show 
that an empirical correction can be used to bring the pro- 
genitor mass distributions in the semi-analytic and N-body 
merger trees into closer agreement. The threshold linear 
overdensity for collapse from the spherical collapse model, 
&c, is replaced with an effective threshold 51^ = fs^dc- In 
the ACDM cosmology employed in the GIF simulation, the 
following form for fg^ was found by Benson et al. to give rea- 
sonable agreement between the progenitor mass functions 
between redshifts and 3: 

fs, = l + O.14[logio(Mhaia//i"'M0) - 15.64], (4) 

where Mhaio is the mass of the final halo at redshift z = 0. 
This form of modification w£is suggested by Tormen (1998). 

These modifications are intended to produce semi- 
analytic merger trees with statistical properties closely 
matched to those of the N-body merger trees. Fig. 4 shows 
the galaxy luminosity functions in the bj and K bands, 
Tully-Fisher relations and global star formation histories for 
both the modified galform model described above (dotted 
lines) and the N-body GALFORM model (dashed lines). It can 
be seen from the figure that these two models predict pop- 
ulations of galaxies with very similar statistical properties. 
The luminosity functions are in reasonable agreement for K 
brighter than about -18 and bj brighter than about -15. The 
Tully-Fisher relations and star formation histories are also 
in close agreement. 

As pointed out previously, the fainter galaxies in these 
models occupy halos with very poorly resolved merger his- 
tories and their properties may be largely determined by the 
effects of limited mass resolution. The solid lines in Fig. 4 
show the properties of the galaxies in the modified GAL- 
FORM model when the minimum halo mass Mres is reduced 
to 5.0x IO^/i-^Mq. This is much less massive than the small- 
est halo Benson et al. were able to resolve in their simula- 
tions and consequently, in this regime, eqn. (4) has not been 
tested and cannot be relied upon to produce a realistic dis- 
tribution of progenitor masses. We do not expect this model 
to reproduce the results of Cole et al. but we show it only 
to provide some indication of the magnitude of the effect of 
introducing low mass halos into the merger trees. 

This "improvement" in mass resolution increases the 
number of faint galaxies, which form in small, previously 
unresolved halos. With a higher minimum halo mass the 
gas in these small halos is unable to cool until it becomes 
incorporated into objects more massive than Mres- This is 
reflected in the luminosity functions which show that there 
are slightly more bright galaxies and far fewer faint galaxies 
at 2: = in the model with poor mass resolution. The star 
formation history is consistent with this, showing that poor 
mass resolution results in reduced star formation at « > 1 
and increased star formation at z « 0. However, calculating 
the global star formation rate involves a sum over all halos. 
At high redshifts this includes a large number of halos of 
low mass whose abundances may be unrealistic due to our 
extrapolation of eqn. (4). Reducing Mres appears to have 
little or no effect on the Tully-Fisher plot. 
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Figure 3. Luminosity functions, star formation histories and TuIIy-Fisher relations for galaxies predicted by the N-body GALFORM 
model using merger trees obtained from the GIF simulation with two different methods of enforcing mass conservation. The solid lines 
show results obtained when mass conservation in the merger trees is enforced by increasing the masses of halos less massive than their 
progenitors. The dotted lines show the results obtained if, instead, the masses of the progenitors of such halos are reduced. 



Overall, the predictions of the N-body GALFORM model 
closely match those of the standard galform model when 
we take into account the differences in the halo mass func- 
tion, the progenitor mass distribution and the mass resolu- 
tion. The differences between the modified GALFORM models 
with high and low mass resolution indicate that, at low lumi- 
nosities, the properties of the galaxies in the N-body model 
are seriously affected by the resolution of the simulation. 
In order to attempt accurately to reproduce the properties 
of observed galaxy populations with bj band magnitudes 
fainter than about -17, an N-body simulation with signifi- 
cantly improved mass resolution would be required. 



4 CONCLUSIONS 

In this paper we have examined how the statistical prop- 
erties of the galaxies predicted by a semi-analytic model 
depend on the way in which the dark matter halo merger 
histories are created. We have developed a method for cal- 
culating merger histories from N-body simulations and used 
the resulting merger trees in a semi-analytic model of galaxy 
formation based on that of Cole et al. (2000). We refer to 
this model as N-body galform and compare it to an oth- 
erwise identical "standard galform " model, which uses 
halo merger histories generated using the Monte-Carlo algo- 
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Figure 4. Luminosity functions, star formation histories and Tully-Fisher relations for three different models. The solid lines correspond 
to the GALFORM model using Monte Carlo generated merger trees as described by Cole et al. (2000), with the modifications explained 
in Section 3 and a mass resolution of 5 x 10^ft.~^MQ. The dotted lines show results from the same model with a mass resolution of 
1.4 X lO^^/i~^M0, equivalent to that of the GIF simulation. The dashed lines show results obtained from the N-body galform model 
which uses merger trees derived from the simulation. 



rithm of Cole et al. This algorithm is bcised on the extended 
Press-Schechter (EPS) theory. 

We find that in a significant number of cases, halos in 
the N-body merger trees are less massive than their progen- 
itors at the previous timestep. When this happens we are 
forced artificially to adjust the masses of the halo or its pro- 
genitors, since in our semi-analytic galaxy formation model 
halos may not lose mass. However, the luminosity function, 
Tully-Fisher relation and global star formation history of 
the galaxies predicted by the semi-analytic model remain 
almost exactly the same whether we add mass to the halo 



or remove mass from the progenitors when we encounter 
this problem. We conclude that the changes we are forced 
to make to the halo masses have very little effect on the 
semi-analytic model. 

If the mass resolution in the standard GALFORM model 
is degraded to that of the N-body simulation and the em- 
pirical fit of Benson et al. (2001a) is used to correct the 
distribution of halo progenitor masses, we obtain luminosity 
functions and Tully-Fisher relations in very good agreement 
with the N-body galform model. This shows that, apart 
from the issue of mass resolution, the only significant statis- 
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tical differences between the N-body merger trees and those 
of Cole et al. are due to the known discrepancy between EPS 
theory and the results of N-body simulations. 

By improving the mass resolution in the standard GAL- 
FORM model to that used by Cole et al. wo were able to 
obtain an indication of the effects of limited mass resolu- 
tion on the N-body model. The mass resolution in the N- 
body merger trees is imposed by the particle mass in the 
GIF simulation, since halos with fewer than 10 particles 
(1.4 X lO^/r-'Mo) arc not resolved. This limitation has a 
noticeable effect on the galajcy luminosity function and we 
find slightly more very bright galaxies, since gas may only 
cool in resolved halos. If only massive halos are resolved, 
cooling is delayed resulting in brighter galaxies at z = 0. 
However, the most obvious effect of poor mass resolution is 
a drastic reduction in the number of galaxies with bj mag- 
nitudes fainter than about -17. This demonstrates that the 
mass resolution of the GIF simulation is insufficient to make 
reliable predictions at these magnitudes. At brighter magni- 
tudes the luminosity functions remain in good agreement. 

In conclusion, when used as the starting point for semi- 
analytic modelling of galaxy formation, merger trees taken 
from an N-body simulation using the technique described in 
this paper result in similar galaxy populations to those ob- 
tained using the (slightly modified) Monte-Carlo algorithm 
of Cole et al. This supports the reliability of our method 
and provides a means to populate large cosmological N- 
body simulations with semi-analytic galaxies at a fraction 
of the computational cost of a hydrodynamic simulation of 
the same volume. When applied to the dark matter com- 
ponent of an Smooth Particle Hydrodynamics (SPH) sim- 
ulation, our model will also allow us to compare SPH and 
semi-analytic treatments of galaxy formation, and in partic- 
ular the cooling of gas within halos, on a halo-by-halo basis. 
This comparison is reported in a companion paper (Helly et 
al. 2002.) 
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